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Top quarks may be produced in large numbers in association with new physics at the LHC. 
The polarization of these top quarks probes the chiral structure of the new physics. We discuss 
several kinematic distributions which are sensitive to the polarization of single top quarks and can 
be used without full event reconstruction. For collimated tops we construct polarization-sensitive 
observables for both hadronic and leptonic decay modes and plot their distributions. We compute 
the observable polarization signals from top quarks produced in the on-shell cascade decay of a stop 
squark into a top quark and a neutralino, as well as top quarks produced in the analogous decay 
chain in same-spin partner models. 



I. INTRODUCTION 

The LHC will provide a wealth of data on the physics 
of the top quark [J Q • Standard Model pair production 
pp — > tt alone is approximately 830 picobarns, allowing 
for detailed study of top properties, such as spin corre- 
lations and rare decays. In addition, any new physics 
which is responsible for electroweak symmetry breaking 
must couple strongly to the top quark, leading to many 
events where top quarks are produced in association with 
new physics. 

Models which aim to stabilize the electroweak hierar- 
chy typically predict a top quark partner, such as the 
stop squark in SUSY or the T" in little Higgs models. In 
order to suppress contributions to precision electroweak 
observables, this top partner is frequently taken to be 
odd under a discrete parity, which in turn typically leads 
to a new stable invisible particle. Thus a standard sig- 
nature of a broad class of well-motivated models is the 
production of a single top quark in the cascade decay of 
an on-shell top partner, accompanied by additional par- 
ticles in the event and possibly missing energy. 

The coupling of top partners to top quarks is in general 
chiral, and a top quark produced in the cascade decay of 
a top partner will be polarized to a degree which de- 
pends on the kinematics of the decay chain as well as 
the couplings. The top undergoes weak decay before the 
process of hadronization can dilute information about its 
polarization, so the angular distributions of its daugh- 
ter particles reflect the polarization of the parent top. 
Observables which measure single top polarization can 
therefore yield measurements which are sensitive to the 
chiral structure of any new physics sector. 

We focus here on kinematic distributions which are 
can be used without fully reconstructed events. To fully 
reconstruct events with missing energy one must make 
assumptions about the event topology and the identities 
of particles not directly observed. This is not always de- 
sirable in events with high multiplicity final states and 



large missing energy. In such events hadronic top decays 
will likely be more useful than leptonic top decays for 
both event selection and top reconstruction, especially 
in the initial stages of analysis. Fully collimated tops 
allow polarization study without full event reconstruc- 
tion in both leptonic and hadronic decay modes, as in the 
collinear limit one may construct observables which be- 
come independent of the unknown boost between the top 
rest frame and the lab frame. 

In section [TT] we discuss several observables whose dis- 
tributions can measure single top polarization, for both 
collimated and uncollimated top quarks. In section IIIII 
we compute the polarization of top quarks produced from 
the cascade decay of a top partner to a top quark and 
an additional particle which is not directly observed. In 
section ITV| we relate this polarization to experimentally 
observable signals, and in section [V] we conclude. 



II. KINEMATIC DISTRIBUTIONS FROM 
STANDARD MODEL TOP DECAY 

In the Standard Model the top decays almost entirely 
through t — > Wb — > fifib- The angular distributions of 
the final state particles depend on the spin of the top. In 
general the angular distribution of any one of the top's 
three daughter particles can be written 



1 oT 1 
TdcosOf ~ 2 ^ 



PtKf COS0f) . 



(1) 



'Electronic address: jshelton@physics.rutgers.edu 



where Of is the angle between the daughter fermion / and 
the top spin axis, as measured in the top rest frame, and 
Vt is the polarization of the top along that axis. We use 
conventions where a top with spin ±1/2 has Vt = ±L 
The "spin analyzing power" of the fermion /, Kf, is a 
calculable coefficient, which depends on the identity of 
/. Unitarity restricts \nf\ < 1. 

The anti-fermion from the W decay — the charged lep- 
ton, in leptonic W decays — is maximally correlated with 
the top spin, kj 2 = 1 [3j. Therefore the charged lep- 
ton is the natural first object to consider in a study of 
top polarization. However, in many inclusive final states, 
pp — ► tX, the leptonic decay mode may be of limited util- 
ity, as the neutrino makes top reconstruction indirect. 
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While it is possible to study top polarization without re- 
constructing the top rest frame, for instance with highly 
boosted tops, the most direct window into top polariza- 
tion uses knowledge of the top rest frame. Moreover, de- 
pending on the particular study, top identification may 
proceed only via hadronic decays where three jets can 
be combined within the appropriate invariant mass win- 
dow. It is therefore interesting to consider polarization- 
sensitive observables that can be constructed for hadronic 
top quarks. 

Hadronically decaying tops present three natural can- 
didate objects whose angular distributions measure the 
polarization of the parent top. One is the 6-jet. At tree 
level, the distribution of the &-quark is given by equation 
(HJ) with 



^m 2 w 



2m 2 w 



-0.4. 



(2) 



One can also consider the angular distribution of the (re- 
constructed) W boson, which is of the form (|TJ) with 



FIG. 1: The distribution 1/YdT/dz of the fraction of visi- 
ble lab frame energy carried by the fe-jet in a highly boosted 
hadronic top, z = £b/£t, plotted as a function of z. The blue 
line (weighted towards high values of z) is the distribution for 
a negative-helicity top quark; the red line (weighted towards 
low values of z) is the distribution for a positive-helicity top 
quark. 



Kw — ~Kb — 0.4. 



(3) 



The spin analyzing power of the 6-jet (or the recon- 
structed W) is less than half of the spin analyzing power 
of the anti-fermion from the decay of the W. Unfortu- 
nately, in hadronic decays, the down-quark jet cannot be 
distinguished from the up-quark jet. However, in the top 
rest frame the down quark is on average less energetic 
than the up quark. Thus the less energetic of the two 
light quark jets in the top rest frame constitutes another 
interesting object whose angular distribution can serve 
to measure the top quark polarization. At tree level, this 
jet has a net spin analyzing power 0, Q 



0.5, 



(4) 



where, again, j denotes the less energetic of the two light 
quark jets in the top rest frame. Leading QCD correc- 
tions to Kb and Kj have been calculated, and are of order 
a few percent Q. In both cases the effect of the QCD 
corrections is to decrease |k|. 

Depending on the method of top identification used in 
a given study, the angular distributions of one or several 
of the three objects mentioned here (the 6-jet, the re- 
constructed W, or the softer of the two light quark jets) 
can form attractive observables which are sensitive to the 
top quark polarization along the top quark direction of 
motion in the lab frame. 



A. Collimated top quarks 

Suitably massive new physics may produce daughter 
top quarks which are highly boosted in the lab frame. 
The top decay products will then be partially or fully 
collimated, necessitating novel approaches to top iden- 
tification and reconstruction [I, 0, i|. When the tops 



are highly collimated, the angular distributions discussed 
above cannot be well measured and other kinematic vari- 
ables are more appropriate. 

We first consider collimated hadronic tops. When the 
6-jet (or equivalently the W) can be separately identified 
within the top, then the natural polarization-sensitive 
observable is the fraction of the lab-frame top energy 
carried by the 6-jet, 



8b 



(5) 



where £/ denotes the energy of a particle / in the lab 
frame. Taking the axis of polarization to be given by the 
top direction of motion, the variable z is simply related 
to cos 9b, 



cos ( 



2m 2 



mf 



(6) 



where (3 is the boost between the top rest frame in the 
lab frame (we work to tree level and in the narrow width 
approximation). The distribution of z is given by 



1 dT 

r dz 



-(1 - K h Vt) + K b V t 



2m 2 



(7) 
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in the collinear limit /3 — > 1. This distribution is plotted 
in Figure [1] for Vt — ±1- Exactly analogous distributions 
pertain for the reconstructed W, interchanging the roles 
of positive and negative helicities. The distribution of 
the total pr of the 6-jet has also been proposed as a 
polarimeter for highly boosted top quarks §]. 

In the highly relativistic limit, leptonic tops can be 
used to study polarization without any need for addi- 
tional event reconstruction, as observables can be con- 
structed which in the collinear limit are independent of 
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the unknown boost between the top rest frame and the 
lab frame. Collimated leptonic top quarks allow the en- 
ergy of the lepton and the &-jet to be separately mea- 
sured. The momentum fraction carried by the lepton 
provides a natural polarimeter. Consider the observable 



the fraction of the visible lab- frame energy carried by the 
lepton. This can be expressed in terms of top rest-frame 
quantities as 



(8) 



2mtE t (l + cos 9 e ) 



2m t Eg(l + pcosQf) + (mf — m^)(l + P cos 6^ cos( -I- p sin^ sin ( cos a) ' 



(9) 



where 0£ is the angle between the lepton momentum and 
the top direction of motion, Eg is the lepton energy in 
the top rest frame, a describes the orientation of the b-£ 
plane relative to the top direction of motion, and 



cos£ 



E e {m 2 t 



E e (m 



(10) 



is the angle between the lepton and the b momenta in 
the b-l plane. The (unknown) boost from the top rest 
frame to the lab frame is parameterized by /?, which we 
will henceforth take to unity. In these coordinates, the 
matrix element for top decay is 



\M\ 2 oc (m t E e - 2E' 2 ){1 +T t cos 6 e ), 



(11) 



where again we take the axis of spin quantization to lie 
along the direction of top motion. The distribution of u 
can then be obtained from 



1 dT 

r du 



dEid cos Qt da \M\ 2 {V t ,E e , cost 
x 6(u — u(Ee, cos 9i, a)). 



(12) 



We use the delta function to carry out the integral over 
a, and evaluate the remaining integrals numerically. The 
resulting distributions for positive- and negative-helicity 
top quarks are plotted in Figure O The distributions 
have a shoulder at 



2m t E, 



2m t Ei 

.rain 

+ (to 



0.215, (13) 



where Ei :Tn i n = j (2mt) is the minimum possible 
value of the lepton energy in the top rest frame. Above 
this value of u, the integrand can receive contributions 
from multiple angular configurations as E( varies, while 
below this value of u, only limited angular configurations 
can contribute. 



III. SINGLE TOP POLARIZATION FROM 
CASCADE DECAYS 

Cascade decays of top partners to tops + missing en- 
ergy are a prototypical example of the processes for which 




FIG. 2: The distribution 1/T dF/du of the fraction of visible 
lab frame energy carried by the lepton in a highly boosted 
leptonic top, u — £e/(£e + £b), plotted as a function of u. The 
blue curve (peaked towards low values of u) is the distribu- 
tion for a negative-helicity top quark; the red curve (weighted 
towards high values of u) is the distribution for a positive- 
helicity top quark. 



the above observables are well-suited. In this section we 
compute the polarization of top quarks produced from 
the decay of a top partner into a top quark and a bo- 
son partner, for both SUSY models (where the boson 
partner is a neutralino) and same-spin partner models 
(where the boson partner is a massive vector boson). We 
work to tree level and in the narrow width approximation, 
and take the top quark to have standard model decays. 
The polarization (Vp) is computed along the axis deter- 
mined by the direction of the top motion in the parent 
rest frame, which is not directly observable; in section HVl 
we will show how to relate this to the polarization (Vd) 
along the axis determined by the direction of the top mo- 
tion in the laboratory frame, which yields predictions for 
experimentally observable distributions. 
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A. Production polarization in SUSY decay chains 

Consider the two-body decay of an on-shell stop squark 
t into a top quark and a neutralino % . The coupling for 
this process can be parameterized as 



att L X +btt R f + H.c, 



where the Yukawa couplings a and b are in general not 
equal. The top quark is thus produced with some nonzero 
polarization, the magnitude of which depends on (1) the 
masses of the particles in the decay chain and (2) the 
mixing of the stop squarks and the neutralinos. We take 
the axis of spin quantization to lie among the top direc- 
tion of motion in the parent stop rest frame. 

The production amplitudes for positive and negative 
helicity top quarks depend on the following functions of 
the top and neutralino energy-momenta, 



F± 



(E t + m t ± \p t \)(E x +m x ± \p x 
y/A(E t + m t ){E x +m x ) 



where all quantities are to be evaluated in the stop rest 
frame. For this on-shell process, the top and neutralino 
energy-momenta are fixed in terms of the masses of the 
particles in the cascade. The functions F± result simply 
from explicit evaluation of the usual spinor wavefunction 
operators y/p ■ a, y/p ■ a on helicity eigenspinors. They 
measure the separate contributions of the two helicity 
states in the left- or right-handed fermion wave function. 
The matrix element to produce a positive (negative)- 
helicity top quark is then the sum of the contributions 
from the left- and right-handed tops, 

M± = Ti(aF T + bF±). 

For finite m t this gives, as expected, a nonvanishing am- 
plitude for top quarks of both hclicities even in the limit 
of a purely chiral vertex, a — ► or b — ► 0. We find the 
net polarization along the production axis to be the 



(\b\ 2 



)M\ Pt \ 



(|o| 2 + \b\ 2 )(ME t - ml) + 2Re(ab)m x m t ' 



(14) 



where M is the mass of the initial squark. The net po- 
larization (Vp) arising from this cascade decay is plotted 
in the solid curves in Figure [3] as a function of the masses 
of the stop squark and the neutralino. For simplicity we 
take purely chiral couplings (6=1, a — 0) in the figure, 
so the plotted polarization is an upper bound. 



B. Production polarization in same-spin partner 
model decay chains 

We now consider the analogous decay chain in the case 
where the top partner T has spin 1/2 and the vector bo- 
son partner has spin 1. We parameterize this process 
by 

-C mt = (at^V L f + bt^V R f)% + H.c. 



-UiU_ — — — SSO 



FIG. 3: Average top quark polarization (Vp) as a function 
of parent particle mass M in GeV, plotted for purely chiral 
couplings (b = 1, a = 0). The red (upper) curves have a fixed 
boson partner mass of 200 GeV. The blue (lower) curves 
have a boson partner mass of M — 200 GeV. Solid (dashed) 
lines plot the average polarization for opposite- (same-) spin 
partners. 



Again, the couplings a and b are in general unequal and 
depend on the mixings of the top and vector boson part- 
ners. The computation is analogous to the one in the 
previous subsection, except now we must also average 
over the spins of the top partner T and sum over the 
spins of the vector boson partner V. We find the net 
degree of polarization along the production axis to be 



(Pp) = 



(\b\ 2 -\a\ 2 )p t (M 2 + 2m 2 v ~m 2 t ) 
(|a| 2 + \b\ 2 ){?>E t m 2 v + 2Mp\) - 6Re(a6)m 2 / m t : 



(15) 

where all kinematic quantities (E t ,p t ) are to be evalu- 
ated in the parent T rest frame. A similar computation 
was carried out in [l(]| for a benchmark model point. The 
net polarization (Vp) arising from this cascade decay is 
plotted in the dashed curves in figure H as a function 
of the masses of the top partner and the vector boson 
partner. The kinematic suppression of the "wrong" po- 
larization state is less for the fermion-fermion-vector in- 
teraction than it is for the fermion-fermion-scalar interac- 
tion. That is, for a given set of masses and couplings, the 
top polarization predicted by a same-spin cascade is less 
than that predicted by an opposite-spin cascade. This is 
due to the additional suppression of the "wrong" helicity 
contribution which comes from the spinor wave function 
of the neutralino in the opposite-spin cascades; in the 
same-spin case, there is some kinematic suppression of 
the "wrong" helicity contributions coming from the lon- 
gitudinal polarizations of the vector boson partner, but 
no suppression from the transverse polarizations. 
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IV. OBSERVABLE TOP POLARIZATION 
SIGNALS 

In events where the top is reconstructed but the rest of 
the event is not, the natural axis for the top polarization 
is the top's direction of motion in the lab frame, which 
we will refer to as the detection axis. The most natural 
experimental observable for hadronic tops is then the an- 
gular distribution of the direction of the bottom jet (or 
the less energetic light quark jet) with respect to the top 
quark direction of motion in the top rest frame, 

The observed polarization is thus along the detection 
axis, (7 7 d). For highly boosted tops, the natural polar- 
ization axis is again the top direction of motion. Notice 
that there is no interference between the two top spin 
states, as the observables we consider are insensitive to 
rotations around the top direction of motion. 

In order to relate the observed polarization to the po- 
larization along the production axis (Vp) computed in 
the previous section, we need to relate the two spin bases. 
The detected polarization (Vd) is related to the produc- 
tion polarization (Vp) by 

(V D ) = (V P ) cos lj, (17) 

where cosw is the angle between the production and de- 
tection axes in the top rest frame. This angle is the 
Wigner angle determined by the composition of boosts 
AtL from the top rest frame to the lab frame, followed 
by Alp from the lab frame to the parent rest frame, and 
finally Apt from the parent rest frame to the top rest 
frame, 

M(A tL )M(A L p)M(A Pt ) = K( u ), (18) 

where A4(A),TZ(9) are the representation matrices for 
the Lorentz transformations [9(. From this composition 
of boosts we find 

E t /3 cos8+ \p t \ 

cosw = — , 

y/ \ Pt | 2 ( 1 + (3 2 cos 2 9)+m 2 f3 2 + 2E t \ Pt \(3 cos 9 

where E t ,pt are the energy and momentum of the top 
in the parent rest frame, 9 is the angle of the top in the 
parent rest frame, and (3 is the boost between the parent 
rest frame and the lab frame. The energy and momentum 
of the top are fixed by the masses of the particles in the 
decay. The phase space integration over cos (9 can be 
carried out simply as (after averaging over top partner 
and vector boson partner spins) the matrix element is 
independent of cos9. This leaves the boost (3 between 
the parent rest frame and the lab frame. In the limit 
where the top partner is produced near threshold with 
little velocity in the lab frame, j3 — > and cos to — > 1, 
as there is little difference between the production and 
detection spin axes. The Wigner angle is also lessened 




FIG. 4: Polarization reduction factor cos a; for three different 
sets of new particle masses, as a function of parent top part- 
ner boost (3. The phase space integral over dcosO has been 
performed. The red (central) curve is for a top partner with 
mass 500 GeV and boson partner 150 GeV. The green (up- 
per) curve is for a top partner with mass 900 GeV and boson 
partner 300 GeV. The blue (lower) curve is for a top partner 
with mass 900 GeV and boson partner 700 GeV. 

when the top quark is highly boosted in the top partner 
rest frame. Thus, as expected, the maximal observable 
polarization signals are realized when the top partner 
is heavy and the vector boson partner is light. We plot 
cos u> as a function of /3 in three different scenarios for the 
new particle masses in Figure [4] (we integrate over cos 9). 
For boosts j3 ~ 0.6, we find that the scenarios with a 
moderate to sizable mass splitting have cosw ~ 0.85 or 
better (the red and green curves shown in the figure), 
while in the case of heavy, degenerate standard model 
partners (the blue curve) the suppression is much larger, 
cosw ~ 0.5. For smaller boosts (3 ~ 0.25 the situation 
is much better, with a suppression cos a; ~ 0.9 even for 
the heavy, degenerate scenario. In principle for a given 
model the distribution of boosts T>((3) can be computed 
from the production dynamics and the PDFs, and the 
observable polarization is then given by 

(V D ) = (Pp)^[ dcos9d[3V(f3) cosv(6,j3), (20) 

allowing detailed quantitative predictions. In practice, 
for sufficiently heavy top partners (and a sufficiently large 
mass splitting between the top partner and its daugh- 
ter particles) , the detailed boost distribution may not be 
necessary. 

V. SUMMARY 

Top quarks produced in association with new physics 
can serve as an important window into the chiral struc- 
ture of physics beyond the Standard Model. We have 
computed here the expected top polarization arising from 
a well-motivated class of events, namely on-shell cascade 
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decays of a top partner into a top quark and a vector bo- 
son partner. Observation of this signal constrains both 
the masses of the particles in the cascade and the chiral- 
ity of the couplings. 

We have focused on variables which can be used in 
events which are not fully reconstructed. In the happy 
circumstance that the LHC reveals a rich spectrum of 
new particles with complicated cascade decays, this strat- 
egy is likely to prove the cleanest way to assemble a large 
sample of inclusive top + new physics events, particularly 
if the new physics includes a stable invisible particle. We 
have discussed several angular distributions relevant for 



hadronic top quarks, and for collimated tops have con- 
structed variables for both hadronic and leptonic decay 
modes and plotted their distributions. 
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